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Abstract 



As follows from the energies of single-particle states in 40 Ca, 90 Zr and 208 Pb nu- 
clei the contribution of many-particle NN forces to the nuclear single-particle poten- 
tial is at least the sum of repulsive and attractive parts resulting from three-particle 
and four-particle forces respectively. In addition the specified nucleon density dis- 
f^*) | tributions in the above nuclei are determined from both the 1 GeV proton-nucleus 

OS ■ elastic scattering and the single-particle energies. 

Os 



1 Introduction 



The nucleon-nucleon interaction proceeds via the meson exchange, and therefore the ex- 
istence of many-particle NN forces is a natural consequence of nonlinearity of the strong 
interaction theory. Indeed, only two-particle forces, Fig. la, are possible in the lowest 
order of linear theory containing only quadratic terms of the meson field if in the meson 
Lagrangian density. The situation is different in nonlinear theories. In the theory with ip 3 
terms the meson may turn into two ones thus giving rise to three-particle forces, Fig. lb, 
in addition to the pairing ones. In the same way the lowest-order four-particle forces, 
Fig.lc, appear in the theory with <f A terms. No higher-power terms exist in renormal- 
izable theories, but the branching of meson is possible in higher orders thus giving rise 
to higher NN forces, Fig. Id. For this reason the empirical information on many-particle 
forces is of fundamental importance for the physics of strong interactions. 

At present the only such information is provided by the calculations for the few-nucleon 
systems. According to these calculations something should be added to the two-particle 
AW forces to get the agreement between the calculated and observed quantities. In 
practice only three-particle forces are added (see Ref . [|l|] and the references therein) since 
the problem is difficult even in this case. Indeed, the forces are characterized by a number 
of adjustable parameters describing the strength, range, the spin-isospin structure etc., 
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and the observed quantities are expressed through these parameters in all orders of the 
perturbation theory. In such conditions it is hardly possible to distinguish, for instance, 
between the genuine three-particle interaction and the combination of the latter with the 
four-particle one. 

As shown in the present paper such possibility may be provided by the experimental 
data on nuclear single-particle states. It is based on the facts that (a) the contribution of 
the two-particle NN forces to the nuclear single-particle potential is fixed by experiment 
and (b) no correlation effects contribute to the single-particle energies [§]. For these rea- 
sons the many-particle contribution may be detected by comparing the observed energies 
of single-particle states with the calculations including the two-particle NN forces only. 



2 Nuclear single-particle states 

We are using the untraditional approach by M.Baranger p|. It is based on the spectral 
representation for the single-particle propagator [||, [|] 

S a (x,x';t) = -i(A o \TiP(x,r)iP + (x',0)\A ) 

; <0(-t) E ^{^t(x'Y- lE]T - tf(r) E Mx)^W)e- lEkT , (1) 



i 
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where \Aq) is the ground-state wave function of nucleus A which is chosen to be even-even 
the ground state thus being nondegenerate, 

^(x) = ((A-l)Mx)\Ao) , E 3 = £,{A)-£ 3 {A-l) 

V fc (a:) = (A \^(x)\(A + l) k ) , E k = £ k (A + 1) - £ (A) (2) 

the sums over j and k thus running over complete sets of states of A — 1 and A+l nuclei. 
The short-time behaviour of a particle or hole which is suddenly created in the ground 
state \A ) is described by the following relations 

EV'jO^/ty) + ^Zi ) k{x)ipt( x ') = i[S A {x,x'; +0) - S A {x,x'; -0) 
j k v 

5(x - x') (3) 
E^'(^)^ + (^) +J2 E kM^k (O = -(s A {x,x'; +0) - S A {x,x'; -0) 

3 k V 

= H sp (x,x') (4) 

J2Efyj(x)ipf{x') +E^(x)^ + (x') = -i(s A {x,x';+0) -S A (x,x';-0) 
j k v 

= H 2 sp (x,x')+U(x,x') , (5) 

where S = dS/dr, S = d 2 S/dr 2 . The single-particle states are treated as eigenstates of 
the single-particle Hamiltonian H sp (x,x'), Eq.(4), 

£\if>\(x) = / H sp (x,x')il) x (x')dx' (6) 
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thus being the doorway states for the single- nucleon transfer reactions As a result of 
the correlation effects these states are distributed over actual states of both the A — 1 and 
A+l nuclei. However the properties of the distribution process permit the determination 
of the single-particle energies from the experimental data. Introducing the spectroscopic 
factors 

2 

(7) 

multiplying Eqs. (3)-(5) by ip^(x)ip\(x r ) and integrating over x and x' we get 

E4 A) +E4 A) = 1 (3a) 

3 k 

E^ A) + E^4 A) = ^ (4a) 

j k 

E^# } + E^4 A) =4 + ^1 (5a) 

j k 

a* = J ^j+( x )U(x,x')ilj x {x')dxdx' . (8) 

The determination of single-particle energies by using the Eqs. (3a)-(5a) and (8) will be 
discussed in Subsect.3.2. 

The Ref. j| derivation of the explicit form of the single-particle Hamiltonian H sp should 
be specified by taking into account the meson-exchange nature of the NN interaction. 

First, the contemporary meson-exchange forces such as the OBE [Q, Paris || and Bonn 
|7j have no hard repulsive core, and therefore the problem of elimination of short-range 
correlations (see the discussion in Ref.fl2|) does not really exist. So there are no actual 
reasons preventing from a direct use of bare NN forces for nuclear structure calculations 
(i.e. the preliminary calculation of the Brueckner G-matrix || or some different way of 
the hard core elimination is not necessary). 

Second, the commutator technique of Ref. |2| does not apply in this case because of the 
retardation. For this reason the field-theoretical approach |3|, |J should be used instead. 
In this approach the single-particle Green function 

G A (x,x';e) = J S a {x,x';t) e l£T dr (9) 

obeys the Dyson equation 

eGa(x, x'; e) = 5(x — x') + % x Ga{x, x'; e) + J M(x,Xi;e)GA(xi,x';e)dx 1 , (10) 

where k x is the kinetic energy and the mass operator M(x, x'\ e) consists of the energy- 
independent part U(x,x'), which is just the nuclear single-particle potential, and the 
energy- dependent one £(x, x';e) which is responsible for the correlation effects: 

M(x,x';e) = U{x, x') + E(x, x'; e) . (11) 



"0a {x)^j,k(x)dx 
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As shown in Refs. j|, |S| the quantity x'; e) vanishes in the e — > oo limit 

E(x,x';£) = (12) 

(the dots in the rhs are the higher-power terms of e^ 1 ) and therefore 

U(x,x') = lim Affox'je) (13) 

e^oo 

the decomposition (11) thus being unambiguous. 

Let us show that the single-particle Hamiltonian of Eq. (4) is 

H sp {x,x') = k x 5(x - x') + U{x, x') , (14) 

whereas the quantity U(x,x') of Eq.(5) is defined by Eq.(12). First let us note that 
according to the time-energy Heisenberg relation the infinite e value is equivalent to the 
infinitely short time interval, the Eqs.(13) and (14) thus meaning that the Hamiltonian 
H sp is indeed responsible for the short-time behaviour of the particle (hole). 

Now let us use the spectral representation |3|, |4[] 

GA(x y., £) = ^^^ + T M^m (1B) 



3 

and the identity 



e — Ej —iS ^ e - E k + i5 



1 1 E E 2 

- + - + — +... (16) 



e-E e(l-E/e) e e 2 e 3 

in the e — >• oo limit. Putting Eq.(16) into Eq.(15) we get the following asymptotic 
expansion 

G A (x,x';e) = h^ + h^ + h^l + ... (17) 

the quantities I (x,x'), Ii(x,x') and I 2 (x,x') being just the lhs of Eqs. (3), (4) and (5) 
respectively. As seen from Eqs. (11), (12) and (14) the Dyson equation (10) may be written 
in the form 

eGa(x, x';e) = 5(x — x') + J ^H sp (x, x\) + _| J Ga(%i, s)dx\ . (18) 

Putting Eq.(17) into Eq.(18) we again get Eqs.(3)-(5), but now the quantities H sp and n 
are determined. 

In terms of bare NN forces (for a moment let us take into account the pairing forces 
only) the single-particle potential is defined by the first-order diagrams of Fig. 2 provided 
the nucleons interact without the retardation 0. This is however not the case for the 
meson-exchange forces including both the momentum and the energy transfer. As a result 
of the latter the exchange diagrams of Fig.2b have the e~ x asymptotics thus contributing 
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to the quantity Ti(x,x';e) rather than U(x,x') [JTIJ. Let us illustrate this for the Bonn B 
potential 0. It is the sum of the terms 



A 2 -fi 2 



2a 



A 2 + q 2 - to 2 



jj? — q 2 — cu 2 



(19) 



in the four-momentum space, the form of the meson- nucleon vertices being specified by 
the Lorentz symmetry of the mesons. The latter is however disregarded because it is 
irrelevant for the energy dependence. With this remark the expression for the Fig.2b 
diagram becomes 



MJx,x'-e) 



d 3 q 



3 iq(r-r') 



iduo 



v(q,uj)GA(x,x';e + uj) 



(2?r) 3 J 2tt 

Considering the case of the monopole formfactor, a = 1, and using the identity 



(20) 



A 2 -/i 2 



A 2 + q 2 — uj 2 J /i 2 + q 2 — uj 2 
d \ 



— lim < ■ 



1 



1 



i i M 



1 _ (A 2 _ ^ 2 ) 

\ dA 2 J 2u A (q) \uj — oJ A (q) + i5 uj + uJ A (q) — i5 / 

M<?) = G" 2 + <7 2 ) 1/2 , Mq) = (A 2 + <? 2 ) 1/2 



(21) 



we get 



M e (x,x';e) = g 2 
-fl-(A 2 -V 



3 iq(r-r') 



d 3 q 

d \ 1 
dA 2 (q) 2uj A (q) 



y e - Ej + u^q) y e-E k - u^q) 



E 



^j(x)^(x') 
e - Ej + u A (q) 



E 



e — E k - uj A (q) 



(22) 



The sign is also irrelevant because the Bonn B potential is the sum of terms with different 
signs. In the e — > oo limit we get 



M e (x,x';e) 



g 2 S(x-x') i 2 
4tt 2 e 1 "' 



_L 1 A 2 -/i 2 

w M (g) wa(?) 2 ^a(<?) 



(23) 



The integral is convergent because the integrand has the g -3 asymptotics. In addition 
to Fig. 2b the quantity Yl(x,x';e) includes an infinite sum of higher-order Feynman dia- 
grams describing all kinds of the correlation effects (Pauli, particle-particle, particle-hole, 
ground-state etc.). 

Thus the only contribution to the nuclear single-particle potential is provided by the 
Hartree diagrams of Fig. 2a and the Hartree-like ones of Fig.3 resulting from the many- 
particle NN forces, the single-particle energies thus being free of the correlation effects. 
The latter is the feature of the Ref.|2| approach, thus permitting the model- independent 
studies of nuclear structure. 
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3 Results 



3.1 Single-particle potential 



The Fig. 2a contribution to the single-particle potential is the convolution of the two- 
particle NN forces with the nucleon density distribution in nucleus. The latter is de- 
termined by the combined analysis of electron-nucleus and 1 GeV proton-nucleus elastic 
scattering data [II] . The two-particle forces are determined from the deuteron properties 



and the elastic NN scattering phase shifts at the energies below the pion production 
threshold. In this way the "pairing" contribution to nuclear single-particle potential is 
fixed by the experimental data thus being independent of any nuclear model. 

We have chosen the Bonn B potential |7|] for the two-particle forces, the motivation 
being discussed in Sect. 4. This choice enables us to check the status of nuclear relativity 
12| by calculating the scalar and vector fields in nuclear interior. Using the parameters 
of Table 5 from Ref.0 and the equilibrium nuclear matter density value p eq = 0.17 fm -3 



we get 10 



V = +284 MeV 



S 



-367 MeV 



These values are close to those provided by the Dirac phenomenology |I3 
tions it is reasonable to treat the single-particle wave functions ip\(x 



(24) 

In such condi- 



as Dirac bispinors 

obeying the Dirac equation. The single-particle Hamiltonian resulting from the Bonn B 
potential is [|Kj 



H sp = -ry°7V + z$(r) 7- + (7 - l)m + V(r) + j°S(r) 



(25) 



m is the free nucleon mass, 7 and 7 are Dirac matrices. The scalar field S(r) consists 
of the isoscalar and isovector parts resulting from the exchange by a and 5 mesons. 
In addition to the isoscalar and isovector parts resulting from the exchange by to and 
p mesons the vector field V(r) includes the Coulomb potential. The small isovector 
quantity $(r) results from the tensor part of the p meson-nucleon coupling. The details 
of the calculations are described in Ref.]l0[|. They also may be found in the extensive 
literature on the Walecka model 



I2H, see for instance Ref. [114" 



The Dirac equation H sp ip\ = e\ip\ is equivalent to Schrodinger-like one for a particle 
with the effective mass 

M(r) = m + -(s{r) -V(r)j (26) 

in the central 

U(r) = V(r) + S(r) (27) 

and spin-orbit potentials [[15[ 

1 d ( 1 



U, 



r dc \2M(r) 



to- . 



(2f 



The terms arising from the quantity <3>(r) are omitted here, but actually they are taken 
into account in the calculations. The quantities in Eqs.(24)-(28) should be supplied by 
the subscript "pair" since they describe the contribution from the two-particle forces. 
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The contribution from the many-particle forces is looked for as the following expansion 
in terms of the nucleon density distribution p(r) 



U m (r) 



he 



a<3p 2 (r) + a 4 p 3 (r) + 



(29) 



the p 2 term results from three-particle forces, Fig. 3a, the p 3 one is of four-particle origin, 
Fig. 3b, etc. (the linear density terms enter the two-particle contribution, Eq.(27)). The 
potential U m (r) is assumed to be equally distributed between the scalar and vector fields, 



S m {r) = V m (r) 



U m (r) 



(30) 



The only motivation for this assumption is the aim to have as little free parameters as 
possible. 



3.2 Single-particle energies 

As mentioned in Sect. 2 the single-particle energies may be determined from experiment 
by using the sum rules (3a)-(5a). The most suitable situation is that for the cases when 
the absolute value of the single-particle energy e\ exceeds the width ex a of the distribution 
region. In these cases all states, over which the single-particle one is distributed, belong 
to the same nucleus, and therefore the sum rules (3a)-(5a) are saturated by only one term 
in the lhs, the first for the hole states and the second for the particle ones. Such situation 
occurs for the peaks in the cross sections of quasielastic knockout reactions (p, 2p) and 
(p, pn) |l6j leading to the hole states with 

|e A | > a ax = 20MeV. (31) 

Indeed, according to Eq.(8) the width a\ depends on the single-particle wave function ipx 
rather than the energy e\, thus being roughly the same for all single-particle states. So 
it is reasonable to identify a with the largest observed value <j max = 20 MeV. The latter 
is the width of the peaks corresponding to the isi/ 2 hole states. For these reasons the 
average energies of the above peaks obeying the Eq.(31) condition may be identified with 
the single-particle energies within the experimental accuracy of 2 4- 3 MeV. In this way we 
demonstrated that the Ref . [Q] approach permits the model- independent determination of 
the single-particle energies. It is worth mentioning that the experimental data for the sp 
energies are independent of those for the two-particle contribution to the sp potential. 

We used the facts that the cross section of the quasielastic knockout reaction leading 
to the fixed nuclear state is proportional to the spectroscopic factor of this state, and the 
absolute values of the s-factors are not necessary when all states, over which the single- 
particle one is distributed, belong to the same nucleus (in this case the relative values 
are sufficient). This is however not the case for weakly bound single-particle states with 
|^a| < Omax- Such states are distributed over actual ones of both A— 1 and A+l nuclei, and 
therefore the s-factors from both the pickup and stripping reactions are necessary. But 
the s-factors are determined with a rather low accuracy because of both experimental and 
theoretical ambiguities, and therefore it is unclear how to pin together the s-factors from 
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the two reactions. For these reasons the energies of weakly bound single-particle states 
are yet unknown. One should also bear in mind that the low-lying states of A ± 1 nuclei 
forming the Fermi-surface of the closed-shell nucleus A are Landau-Migdal quasiparticles 
0] rather than the single-particle states of nucleon. Indeed, the correlation term x f ; e) 
is included in the quasiparticle energies in contrast to the single-particle ones, see Ref.[(|] 
for details. 



3.3 Many-particle forces and specified density distributions 

The observed energies of neutron and proton single-particle states in 90 Zr are plotted in 
Fig. 4 together with the results of the calculations. As seen from the figure the lowest 
"pair" states lsi/2 are significantly underbound whereas the higher states, especially the 
2s\/2 ones, are overbound. The same compression of "pair" single-particle states occurs 
in 208 Pb, Fig. 5, where this effect is even more pronounced, and also in 40 Ca, Fig. 6. This 
means that the potential well, including two-particle forces only, is somewhat too wide 
but insufficiently deep. The isoscalar potentials in 90 Zr, "pair", actual and the many- 
particle contribution, are plotted in Fig. 7. As seen from the figure, the many-particle 
contribution (i.e. the difference between the actual and "pair" wells) consists of the 
repulsive and attractive parts, the radius of the latter being less than that of the former. 
So the expansion (29) has to contain at least two terms of different sign obeying the 
above condition. The most simple possibility is provided by the sum of the three-particle 
repulsion, a 3 > 0, and the four-particle attraction, < 0. Of course we cannot guarantee 
the absence of contributions from higher many-particle forces. It only should be mentioned 
that the above possibility corresponds to the least number of free parameters. 

Taking into account the possible contribution of many-particle forces to the isovector 
nuclear potential the quantity U m (r) is chosen as 



U m (r) = hc^a 3 p (r) + a A p (r) - r 3 |a 3 p(r) + a 4 p (r)J p (r) j , 
p(r) = p n {r)+p p (r), p'{r) = p n (r) - p p (r) . (32) 

p n (r) and p p {r) are neutron and proton density distributions in nucleus. As seen from 
Fig. 4 the description of single-particle energies is improved by including the many-particle 
contribution. 

The results, which are labelled as "tot" in Fig. 4, are obtained using the Woods-Saxon- 



like density distributions of Ref. fllTfl which are folded with the nucleon electromagnetic 
form factor, see Ref.[|ll[] for details. However the nuclear potential is expressed through 
the point densities since the finite size of nucleon is taken into account in the NN forces. 
For this reason we used the point densities PqWa{t) which are obtained from those of 



Ref. JTT| by usual deconvolution procedure. 



It should be mentioned that the electron and proton elastic scattering data under- 
lying these densities are sensitive to nucleon density distributions in the surface region 
of nucleus, whereas the single-particle energies are sensitive to those in nuclear interior. 
Therefore the observed single-particle energies may be used to specify the nucleon density 
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Table 1: Neutron and proton density parameters 









a p 


Pp 


90 Zr 


-0.0836 


0.4823 


-0.0287 


0.4416 


40 Ca 


-0.0513 


0.5469 


-0.0215 


0.4557 


208p b 


-0.2146 


0.4670 


-0.0092 


0.2200 




-0.3083 


0.4646 


-0.1308 


0.2779 



distributions. We tried many different forms of p(r). The most appropriate one is found 
to be 

p(r) = Po(w A (r) + aW A (0)MM) , (33) 

where y?4(x) is the fourth order Hermite function. The best fit parameters a and (3 of 
neutron and proton density distributions are shown in Table 1. The best fit strength 
parameters are found to be the same for all nuclei. They are 

a 3 = 13.6608 fm 5 , a A = -80.2568 fm 8 (34) 
a 3 = 28.7531 fm 5 , a7 = -144.3534 fm 8 . 

The calculations with the specified densities are labelled as "totl" in Figs. 4-6. 

Slightly better agreement for 208 Pb is provided by the following strength parameters 

a 3 = 15.1120 fm 5 , a 4 = -90.9870 fm 8 , 

a7 = 21.0000 fm 5 , aj = -99.0890 fm 8 (35) 
with the following form of the nucleon density distributions 

p(r) = Po W A (r)(l + aMPr)) . (36) 

The corresponding parameters a and (3 are those of the fourth row in Table 1, the results 
are labelled as "tot2" in Fig.5. 

The largest discrepancy between the observed single-particle energies and the "totl" 
results for 90 Zr and 40 Ca as well as the "tot2" ones for 208 Pb is less than 3 MeV, the 
average discrepancy is 1.65 MeV. 

The specified nucleon density distributions are plotted in Figs. 8 and 9. As seen from 
the figures the neutron densities have a pronounced dip in the center of nucleus which 
increases with increasing mass number. 

We also calculated the 1 GeV proton elastic scattering cross sections within the 
Glauber-Sitenko theory [[3j]] using both the specified densities and those of Ref . ]ll . The 
results, see Fig. 10, clearly show that the agreement with experiment is equally good for 
both the Ref. JTl]] and the specified densities. 
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The possible reasons for the difference between the Eq.(34) and the Eq.(35) strength 
parameters are as follows, (i) Expression (32) corresponds to the zero-range forces 

f3(r 12 ,r 13 ) = (a 3 + T 1 T 2 a 3 ~)5(r 12 )5(r r3 ) (37) 
fa(ri2,ri3,ru) = (a 4 + t x t 2 o± )5(r 12 )5(ri 3 )5(r 1A ) , (38) 

whereas the actual forces may be of finite range, (ii) The above forces should have been 
folded with the two-particle density and the three-particle one rather than the products 
of single-particle densities, (iii) As mentioned above, the contributions from higher many- 
particle forces may be present. 

The difference is however small for the charge-symmetric nuclear matter with the 
equilibrium density p eg = 0.17 fm -3 . Indeed, as follows from Eqs. (32), (34) and (35) the 
contributions of three (four)-particle forces to the isoscalar single-particle potentials are 

U 3 = hca 3 p 2 eq = 78 (86) MeV , (39) 
U 4 = hca 4 pl q = -78 (-88) MeV . (40) 

The first values in the rhs correspond to the Eq.(34) parameters whereas those in paren- 
theses refer to the Eq.(35) set. They are not small compared to the "pair" value t/ pa i r = 
—83 MeV, see Eqs. (24) and (27). But they nearly compensate each other thus giving 
rise to the conclusion that the isoscalar part of the nuclear single-particle potential is 
mainly of the two-particle origin. This conclusion is supported by the fact that according 
to Eq.(30) 

S 3 = V 3 = 39 (43) MeV , S 4 = V 4 = -39 (-44) MeV , (41) 

thus being considerably less than the pair values S pajir = —367 MeV, Vp air = 284 MeV, see 
Eq.(24). 

The situation is different for the isovector part of the potential. With p~ = ^-j^-p eq 
Eq.(32) gives 

U~ = hc{a^p eq + a- A p 2 eq )p- = 24^— MeV (42) 

for both Eq.(34) and Eq.(35) parameters whereas the "pair" value provided by the Bonn B 
potential is U~ a[l = 7^j^MeV, the isovector nuclear potential thus being mainly of many- 
particle origin. The reason is due to the fact that the two-particle contribution arises 
from the exchange by isovector mesons p and 5 which are weakly coupled to nucleon |7J . 

4 Summary 

In this way we demonstrated that the many-particle NN interaction includes at least 
the three-particle repulsion and the four-particle attraction. This result is restricted 
because no information is obtained for the many-particle forces which do not contribute 
to the nuclear single-particle potential. But it provides a very instructive example of 
the situation where the experimental data on complex nuclei are more appropriate for 
the fundamental problem than those on few-nucleon systems. Indeed, the latter ones 
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are described by solving the complicated many-particle quantum mechanical problem 
including the interaction in all orders of the perturbation theory. On the contrary the 
single-particle states are solutions of the one-particle problem which is much more simple. 
In addition the nuclear single-particle potential is expressed through the NN interaction 
in first order of the perturbation theory, the result thus being visual (see Fig. 7). 

Besides the above-mentioned A dependence of the strength parameters our results 
may have the following additional ambiguities. 

1. They are essentially based on the choice of the Bonn potential for the two-particle 
NN forces. The reason for this choice is a very high level of confidence: (a) the physics 
underlying the Bonn potential is absolutely clear, (b) it contains only one adjustable 
parameter, (c) the two-nucleon data are described with x 2 /datum=1.9 [0 (according to 
this reference the Bonn B version we used in the present work is completely equivalent 
to the full one). 

Of course our results may be changed with the progress of the knowledge about the 
NN interaction. We insist however that the acceptable new-fashioned potential must be of 
higher level of confidence than the Bonn one. This means that (a) the underlying physics 
should be as clear as that for the Bonn potential, (b) the number of free parameters may 
be also only one but the x 2 /datum value should be considerably less or (c) the x 2 /datum 
value may be the same but no free parameters should be present. 

2. In the Ref . [|IT)f| experiments the energy of the knocked-out nucleon is only about 
100 MeV. This may be insufficient to neglect the final-state inelastic interactions leading 
to additional excitation of the final nucleus. As a result of such excitations the average 
energies of the peaks may be shifted from the single-particle energy values because the 
reaction mechanism is not a pure quasielastic knockout in this case. To get rid of this 
ambiguity the additional quasielastic knockout (p, p'N) or (e, e'N) experiments are desired 
in which the energy of the knocked-out nucleon would be about 1 GeV. We hope that our 
work will stimulate such experiments. 

The authors are indebted to Dr. M.B. Zhalov for valuable discussions. 
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Figure Captions 



Fig.l. Two-particle (a), three-particle (b), four-particle (c) and higher (d) nucleon-nucleon 
interaction forces. The wavy lines are for mesons, the full ones are for nucleons. 

Fig. 2. First-order Hartree (a) and exchange (b) diagrams. 

Fig. 3. Contributions of three-particle (a) and four-particle (b) forces to nuclear single- 
particle potential. 

Fig. 4. Single-particle energies of neutron and proton states in 90 Zr. The energy scale is 
shown from the left of each figure. The labels are " exp " for the observed energies, 
"pair" for the calculations taking into account the two-particle forces only, "tot" for 
those including the many-particle forces and using the density distributions from 
Ref.|lll], and "tot 1" for the case when both the many-particle forces and specified 
densities are included. 

Fig. 5. The same for 208 Pb. The label "tot 2" is for the case when the individual strength 
parameters for 208 Pb are used. 

Fig. 6. The same for 40 Ca. 

Fig. 7. Isoscalar potential in 90 Zr. The full, dashed and dot-dashed lines are for actual, 
"pair" and many-particle part of the potential respectively. 

Fig. 8. Neutron density distributions in 40 Ca, 90 Zr and 208 Pb. 

Fig. 9. The same for protons. 

Fig. 10. Elastic scattering of 1 GeV protons on 40 Ca, 90 Zr and 208 Pb. The calculations 
with the specified densities are shown by full line, those with the Woods-Saxon-like 
ones from Ref. [11] are plotted by dashed line, the dots are for the experimental data. 
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